Spin density in frustrated magnets under mechanical stress: Mn-based antiperovskites 
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In this paper we present results of our calculations of the non-colUnear spin density distribution 
in the systems with frustrated triangular magnetic structure (Mn-based antiperovskite compounds, 
MnsAN (A=Ga, Zn)) in the ground state and under external mechanical strain. We show that 
the spin density in the (lll)-plane of the unit cell forms a "domain" structure around each atomic 
site but it has a more complex structure than the uniform distribution of the rigid spin model, i.e. 
Mn atoms in the (lll)-plane form non-uniform "spin clouds", with the shape and size of these 
"domains" being function of strain. We show that both magnitude and direction of the spin density 
change under compressive and tensile strains, and the orientation of "spin domains" correlates with 
the reversal of the strain, i.e. switching compressive to tensile strain (and vice versa) results in 
"reversal" of the domains. We present analysis for the intra-atomic spin-exchange interaction and 
the way it affects the spin density distribution. In particular, we show that the spin density inside 
the atomic sphere in the system under mechanical stress depends on the degree of localization of 
electronic states. 

PACS numbers: 



I. 1. INTRODUCTION 

In magnetic materials the local magnetic moments 
(LMM) are generally due to the d- and f- electrons. These 
electrons are localized, for example for 3d-metals the ra- 
dius of d-orbitals is less than lk{localized region (LR)). 
On the other hand, s- and p-orbitals in 3d-metals are less 
localized or even delocalized {delocalized region (DR)). 
Thus, there are regions of high spin density (SD) in the 
LR near the nucleus and low SD further away in the DR. 
Despite the inhomogeneity of SD in magnetic materials, 
it is frequently described within the quasispin approx- 
imation [1], according to which the direction of the SD 
around each atom is taken constant within atomic sphere 
(AS) or polyhedron. This approach works better for the 
systems with strong spin-exchange coupling because the 
strong exchange between spins ensures that the SD di- 
rection stays almost uniform. Yet, for some systems it is 
essential to consider detailed distribution of magnetiza- 
tion axis as this may reveal additional important features 
of the magnetic configurations [2, 3] overlooked within 
the rigid spin approximation (RSA). In this paper we ex- 
amine the SD distribution of the systems with frustrated 
magnetization, in particular MnaAN antiperovskites, us- 
ing fully unconstrained SD functional without any as- 
sumption on the uniformity of SD around each atom. 

Many interesting phenomena in magnetic materials are 
due to the interplay of magnetic and structural degrees 
of freedom because the SD is sensitive to mechanical de- 
formations such as compression, in-plane biaxial strain, 
etc. In particular, such phenomena based on magneto- 
mechanical coupling have been reported for the Mn-based 
antiperovskites, - for example piezomagnetic effect [4], 
a giant magnetoresistance of more than 10% in pulsed 
magnetic fields [5] , invar effect (or even negative thermal 
expansion) [6]. While the SD evolves under mechanical 
stress, the RSA assumes that its direction is changing 
uniformly inside of the AS. Yet, it is not obvious that in 



the excited states the spin moments in LR and DR rotate 
to the same angle. We demonstrate in this paper that the 
spin directions in these two region may rotate in opposite 
directions. 

The ground state of MnsAN is shown in Fig. [1] (top 
panel). It forms an antiperovskite crystal structure with 
cubic space group PmSm. The orientation of the LMMs 
of the Mn atoms forms a non-collinear, F^^ structure 
in the Bertaut's classification [7]. This is a structure 
with the spins on the (lll)-plane with clockwise or coun- 
terclockwise spin configurations, such that the spin mo- 
ments in the plane are canceling each other. This some- 
what puzzling magnetic structure raises a natural ques- 
tion on how one treats the SD distribution in the system. 
In our previous work [4] we analyzed the piezomagnetic 
properties of MnaAN without presenting any details on 
the SD distribution, but rather we have discussed LMMs 
integrated over the atomic spheres. 

The paper is organized as follows: in Section 2, we 
present results of our calculations of the magnetic prop- 
erties of antiperovskite such as MnaZnN under external 
stress; in Section 3, we present our results on the SD 
distribution in the (lll)-plane of the unit cell, and we 
show that it has more complex structure than the one 
postulated in the RSA; in Section 4 we analyze the struc- 
ture of the SD distribution in terms of the intra-atomic 
spin-exchange interaction; and finally, the last section is 
devoted to the conclusions. 

We employ the projector augmented- wave (PAW) 
method originally proposed by P. Blochl Q. We use 
its implementation by G. Kresse and D. Joubert in the 
Vienna ah initio simulation package (VASP) code 
within a Perdew-Burke-Ernzerhof (PBE) generalized gra- 
dient approximation [TO] of the density functional theory 
(DFT). We use al2xl2xl2 k-point sampling and the 
Blochl's tetrahedron integration method [11]. We set the 
plane-wave cut-off energy to 300 eV and we choose the 
convergence criteria for energy of 10~^ eV. Within this 
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FIG. 1: (color online) Top panel: non-collinear F^^ structure, 
unit cell. Blue spheres - Ga atoms; red spheres - Mn atoms; 
green sphere - N atom. Arrows represent LMMs of Mn atoms. 
Bottom panel: the net magnetic moment per unit cell and the 
rotation angle of Mn2 LMMs as functions of biaxial strain 
(Aa/a (%)) in MnsZnN. Inset on bottom panel: schematic 
view of the variation of Mn LMMs for MnsAN in (lll)-plane 
as a function of strain. Yellow arrows - compressive strain, 
black arrows - tensile strain, red arrows - ground state. 



framework we perform the non-collinear magnetic struc- 
ture calculations, and we allow the LMMs to relax to the 
equilibrium configuration. 



II. 2. MAGNETIC PROPERTIES OF 
MN-BASED ANTIPEROVSKITES UNDER 
STRESS 

The equilibrium parameters for the MnaZnN are as 
follows: 3.87A for the lattice constant and 2.5 /iB for 
the LMM integrated over the volume of the space-filling 
AS. The magnetic moment is extremely sensitive to the 
volume change. It steadily increases from 0.5 /j^b to 3.7 
/i^ as the volume increases from the lattice parameter 
3.5A to 4.3A. Application of the external strain results 
in induced magnetization. As one can see from the inset 
on the bottom panel of the Fig. [1] the magnetic moment 
of the Mnl atom in the basal plane (formed by atoms of 
Mn and Ga) preserves its direction but its magnitude de- 
creases/increses under compressive/tensile strain. LMMs 



of the Mn2 atoms in the Mn2N plane change both their 
direction and magnitude, - they rotate towards the [110] 
direction when the compressive strain is applied in such 
a way that they become almost opposite to the LMM of 
Mnl atom at large strains. The LMMs of Mn2 atoms 
become more aligned with each other (see Fig. [TJ - in- 
set on the bottom panel). When tensile strain is applied 
the LMMs of Mn2 atoms rotate in direction opposite to 
the one at the compressive strain and the Mn2 moments 
become more anti-aligned. The magnitude of the LMMs 
generally decreases with the compressive strain and in- 
creases with the tensile strain. Thus, MnaAN (A = Ga, 
Zn) antiperovskite acquires a net magnetization under 
biaxial strain as a result of the rotation of LMMs and 
change in their magnitude. The net magnetization is 
directed along the C22 axis (this corresponds to a diag- 
onal axis C2a in the cubic cell). Fig. [1] (bottom panel) 
shows magnetization per cell (black line with solid cir- 
cular data points) as a function of biaxial strain. The 
strain induced magnetization is about 0.07 /j^b per unit 
cell for the MnsZnN at 3% of the compressive strain (for 
MnaGaN this value is substantially larger, in particu- 
lar it is about 0.16 fiB at 2% of the compressive strain). 
The piezomagnetic effect in MnsAN is /mear and exhibits 
magnetization reversal with the applied strain, i.e. the 
direction of the magnetization along the C22 axis reverses 
with the reversal of the applied strain. The rotation an- 
gle of the LMM of Mn2 atom under biaxial strain with 
respect to the Mn moment direction in the ground state 
as a function of strain is shown in Fig. [1] (bottom panel). 
The dependence is linear in the range of ~ ±2% of strain, 
and the LMMs rotate in opposite direction if the sign of 
the strain is reversed. 



III. 3. SPIN DENSITY DISTRIBUTION 

Fig. [2] shows the SD distribution in the (lll)-plane 
of the unit cell for the non-collinear ground state F^^ 
structure of MnaAN. The overall structure of the SD in 
the (lll)-plane resembles vortex with the center in the 
middle of the triangle formed by the Mn atoms. The 
orientation of the spin density in the center of the vor- 
tex is shown schematically on the inset of the Fig. [2l It 
reflects the anti-ferromagnetic nature of the inter-atomic 
exchange coupling of the nearest neighbor Mn atoms in 
the (lll)-plane. Red arcs on the plot schematically show 
the AS within which the SD direction is assumed to be 
uniform by the quasispin approximation. As one can 
see from the plot, this uniformity is only very approxi- 
mate as the magnetization axis direction clearly changes 
within the "quasispin domain". Particularly, in the LR 
(approximately inside r<0.7A sphere) spin axis direction 
is more or less uniform. However, in the DR of low SD 
(approximately 0.7A<r<1.4A) the direction is changing 
substantially (by about 30 degrees). 

Even more dramatic effect can be observed un- 
der strain, where regions of SD inside AS rotate 
in opposite directions. To better visualize the non- 
uniformity of SD we plot the vector product of the 
magnetizations in the ground state and under strain 
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(m (equilibrium) xm (strain)) in the (lll)-plane (see Fig- 
ure [3]). The out of plane arrows (green) show counter- 
clockwise rotations, while into the plane arrows (blue) 
show clockwise rotations of SD under strain. One can 
clearly visualize from Figure [3] the appearance of do- 
mains with different direction of SD rotation. Partic- 
ularly, there are relatively large domains in DR which 
rotate opposite to the domains in LR. Therefore, the SD 
distribution in the (lll)-plane is more complex than the 
ideal picture of the uniform SD around each atom as- 
sumed by the RSA. To understand this "beyond RSA" 
pattern we analyze the intra-atomic spin-exchange inter- 
action and the way it affects the SD around each atom. 




Mnl 



FIG. 2: (color online) Spin density distribution in (lll)-plane 
for the non-collinear F^^ structure of MnsAN, unit cell. Red 
spheres represent Mn atoms. Red arcs schematically high- 
light SD distribution "domains" around each Mn atom. Blue 
arrows represent SD direction. Inset schematically shows the 
orientation of the spin density in the center of the vortex. 



IV. 4. DISCUSSION 

To simplify our analysis we consider SD along the 
line between Mn2 and Mn3 atoms, i.e. we reduce 3- 
dimensional problem to 1-dimensional. Fig. [H (top panel) 
shows the ground state SD distribution for the MnaGaN 
along the lines connecting Mn atoms in the (lll)-plane 
of the unit cell. The angle, 9 between the SD direction 
(blue arrows on the plot) and the z-axis (normal to the 
line between Mn2 and Mn3 atoms) does not show the uni- 
formity inside AS around each Mn atom assumed in the 
RSA. Fig. m (bottom panel) shows the magnitude of the 
SD, |m| as a function of distance between Mn2 and Mn3 
atoms. Like ^, the |m| also shows distinct non-uniform 
behavior. Qualitatively these features of SD distribution 
can be interpreted in terms of the spin-exchange interac- 
tions. For continuous SD distribution the spin-exchange 
contribution to the interaction energy between Mn2 and 
Mn3 atoms can be written in terms of the spin (Heisen- 
berg) Hamiltonian, H{x^x') [13] as follows: 




FIG. 3: (color online) Vector product of the local SD in the 
ground state and under strain in the (lll)-plane: MnsGaN 
(a, b) and MnaZnN (c, d): red spheres represent Mn atoms; 
green (light) domains represent out-of-plane vector product 
direction; blue (dark) domains represent into-the-plane vector 
product direction. 



EexX^) = - dxdx' J{x,x')SxSx'Cos{Ox - Ox') =^ 

Jx2 ^ V ' 

H{x,x') 

1 + ' 

E ex. ^ const -\- - J dxdx' J{x^x')sxSx'e ~ {60)'^ 

1 f^^ 

= const ^- / dxJ{x)sle~^{SOf{l) 

2 Jx2 

where J{x^x') is the exchange parameter; Sx and Sx' are 
spin vectors at x and x'] x and x' are "spin positions" 
along Mn2^Mn3 line; x2 and x3 are coordinates of the 
Mn2 and Mn3 atoms; Ox and Ox' are angles between spin 
at X and x' and z-axis. Besides, here we used the ex- 
ponential behavior of the magnitude of the spin vector 
(see Fig. m bottom panel), and J-function-like behavior 
of J(x, x')^ i.e. 

SxSx' — SxSx'C ^ (2) 



J{x^x')sx'dx' = J{x)sx (3) 

We also assume that x and x' are very close to each 
other, therefore x ^ x' ^ 2x. 

To find the expression for the 0{x) which fits the SD 
distribution shown on Fig. [H (top panel) we have to min- 
imize Eq(T] w.r.t. X, i.e. we have to set the functional 
derivative of Eqil] equal to zero. This results in the 
following differential equation with the boundary con- 
ditions: 
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FIG. 4: (color online) Magnetization distribution between Mn 
atoms in (lll)-plane for the non-collinear F^^ structure of 
MnsGaN (top panel); |Spin density | as a function of distance 
between Mn2 and Mn3 (bottom panel) for ground state (red), 
1.55% tensile strain (blue), and 1.55% compressive strain 
(black). 



J{x) ' 0{x) + d{x) ■ j{x) = 0; 
0{x2) = 30°; 0{x3) = -30° (4) 

from which 

(P . d ^, . d 
e V2 . -—e{x) + — i9(x) • —e V2 = O; 
a^x ax ax 

0{x2) = 30°; 0{x3) = -30° (5) 

We solve this equation separately for LR and DR by 
choosing the value for A which corresponds to the best fit 
for the distribution of ^ as a function of distance between 
Mn2 and Mn3 atoms. Results of the solutions for the 
ground state and strained MnsGaN are summarized on 
Fig. El (top panel) and Table H 

Fig. [5] (top panel) shows the change in for MnaGaN 
along the half-distance of the line connecting Mn2 and 
Mn3 atoms for the ground state, tensile strain (~1.55%), 
and compressive strain (~1.55%) (for the other half of 
the Mn2-Mn3 distance we just have to change the sign 
of 0). As one can see the angle under tensile strain for 
the SD distribution in the LR is lower than the for the 
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x: Mn2 - Mn3 half-distance (A) 

FIG. 5: (color online) Top panel: ^ as a function of half- 
distance between Mn2 and Mn3 atoms for ground state (red 
spheres), 1.55% tensile strain (blue stars), and 1.55% com- 
pressive strain (black squares). Bottom panel: ^ as a function 
of half- distance between Mn2 and Mn3 atoms. 



J{x) and 


Equilibrium 


1.55% tens. 


1.55% compr. 


J(x)(LR) 


^-x/0.17 


^-x/0.18 


^-x/0.16 


J(x)(DR) 


^-x/OAO 


^-x/0.37 


^-x/0.50 


6'(LR) 


31-0.35-e^-^" 


29-0.40-e^-^"^ 


33-0. 29-e^-^"^ 


0{DR) 


21-0.79-e2-^"^ 


20-0.57-e^-^"^ 


23-1.69-e^"^ 



TABLE L J(x) and 0{x) for the MnsGaN in the ground state 
and under strain for regions of localized (LR) and delocalized 
(DR) electronic states. 



ground state in the same region. This is consistent with 
the general description of the piezomagnetic effect pre- 
sented above, - see for example the inset on the bottom 
panel of the Fig. [H Yet, in the DR the ground state 6 
becomes smaller, right up to the point when changes 
the sign. Under compressive strain the is larger for 
the SD distribution in the LR than the for the ground 
state in the same region (again, consistent with the de- 
scription of the piezomagnetism presented above). Closer 
to the DR the difference in angles for the system under 
compressive strain and for the system in equilibrium de- 
creases and finally changes its sign, after which the for 
the ground state becomes larger than the for the com- 
pressed system. But for both signs of the strain (tensile 
and compressive) the magnitude of the difference in the 

for the strained and equilibrium states (^ =^ (under 
strain) - (equilibrium)\o>oore<o) is larger in the LR com- 
pared to the DR (see Fig. [5] (bottom panel)). This ^ sign 
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changing mechanism is obviously in agreement with the 
SD distribution picture shown on the Fig. [3] (see for ex- 
ample (a) and (b), - both show the reversal of the sign 
of the m(equil.)xm(strain) at the border of the LR/DR 
domains. These additional spin domains'' in DR can- 
not be revealed within the RSA framework. The sign of ^ 
is changing because the strain affects the overlap of the 
wave functions (spin-exchange interaction). For exam- 
ple, when tensile strain is applied the distance between 
Mn2 and Mn3 atoms increases which results in increased 
localization of the electronic states in the LR (in J{x) 
A increases from 0.17 to 0.18, - see Table H]). At the 
same time, J{x) of delocalized states in the DR becomes 
weaker (A decreases from 0.40 to 0.37). This decrease in 
J{x) in the DR results in reversal of the ^ sign, which in 
its turn leads to appearance of the additional " domains" 
on Fig. [3l On the other hand, when compressive strain 
is applied, the states in the LR become less localized 
(A decreases from 0.17 to 0.16) while the J{x) in DR be- 
comes stronger (A increases from 0.40 to 0.50). This J{x) 
increase in DR under compressive strain also results in 
reversal of the ^ sign (as in the case of the tensile strain 
the spins in LR and DR turn in different directions). All 
the above arguments based on our model present gen- 
eral idea of the mechanisms responsible for the "beyond 
RSA" picture. 

V. 6. CONCLUSION 

We have presented results of our calculations on the 
spin density distribution of the MnaAN in the ground 



state and under external bi-axial strain. In both cases 
the spin density demonstrates non-uniform features such 
as appearance of additional spin domains in the (111)- 
plane not revealed within the quasispin approximation. 
The orientation of the spin density depends on the di- 
rection of the strain. We have presented explanation of 
the non-uniform distribution of the spin density in terms 
of the intra-atomic spin-exchange interaction. In par- 
ticular we have shown that the spin density distribution 
depends on the degree of localization of the electronic or- 
bitals, i.e. strongly and weakly localized orbitals within 
the same atomic sphere rotate in different directions un- 
der external strain. 
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